Pair breaking by impurities in the two-dimensional t — J model 
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Pair breaking mechanisms by impurities are investigated in the two-dimensional t — J model by exact 
diagonalization techniques. Analysis of binding energies, pairing correlations, dynamical spin and 
pair susceptibilities shows that non-magnetic impurities are more effective in suppressing pairing 
than magnetic ones in agreement with experimental studies of Zn- and Ni- substituted High-Tc 
superconductors. 
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Many recent experiments of hole-doped cuprate super- 
conductors showing a large anisotropy consistent with 
gap nodes on the Fermi surface jj] have been interpreted 
in terms of an unconventional d-wave pairing state. This 
possibility of d-wave superconductivity has in turn re- 
newed the interest in studying the effects of impurities 
on the superconducting state. It is well-known that mag- 
netic impurities are strong pair breakers for singlet su- 
perconductors while nonmagnetic impurities have a pair 
breaking effect only for higher orbital momentum states 
such as a d-wave pairing state. |2| In fact, in several re- 
cent studies it was suggested the possibility of using the 
response of these materials to different types of impu- 
rities to distinguish between an unconventional e?-wave 
pairing and alternative scenarios such as a more con- 
ventional anisotropic s-wave pairing. |3|-||] Experimental 
studies have shown that divalent Zn and Ni ions go to the 
planar Cu(2) sites with the Zn impurity having a 5 = 
configuration and the Ni impurity a 5 = 1 configuration. 
H The important experimental fact is that the nom- 
inally nonmagnetic impurity Zn is more effective than 
the magnetic impurity Ni in destroying superconductiv- 
ity. In conventional superconductors, with the usual 
s-wave superconductivity, a magnetic impurity is most 
destructive due to the magnetic pair-breaking effect. ||] 
Although the observed behavior could be characteristic 
of a d-wave pairing, it has been also suggested that a Zn 
impurity may not simply behave as a vacancy because it 
could induce a magnetic moment in the CuC>2 planes due 
to the strong correlations present. [pPJlOtl However, other 
experimental studies |[l| have indicated that the esti- 
mated impurity moment- carrier exchange in Zn doped 
YBa 2 Cu307 is too small to account for the Tc suppres- 
sion by a magnetic pair-breaking mechanism. 

From a theoretical point of view it is of extreme impor- 
tance to explain at least qualitatively the effects of both 
magnetic and non-magnetic impurities on the supercon- 
ducting properties of the cuprates as seen experimentally. 
These effects have not been clearly explained by existing 
analytical macroscopical theories. |p|-p|,^2|-p^| An alter- 
native approach is to analyze them in the context of mi- 



croscopic models using numerical techniques. This pro- 
gram has been recently initiated by Poilblanc, Scalapino 
and Hanke |l(]Jl7|] by studying the two-dimensional (2D) 
t — J model in the presence of a single impurity. The 
two-dimensional t—J model has been extensively stud- 
ied in the context of high-T c superconductivity since it 
contains the essential low-energy physics of the CuC>2 
planes present in the cuprates. |L8| The 5 = Zn impu- 
rity was modelled by an inert site, |16| while the 5 = 1 
Ni impurity was approximated by a static spin-1/2 (see 
Hamiltonian (Q) below). Jl|] The main results obtained 
by these authors is that a hole binds to the impurity with 
different spatial symmetries. 

The purpose of this Letter is to study the effect of 
magnetic and nonmagnetic impurities on several quanti- 
ties related to pairing of holes in the 2D t—J model near 
half-filling. This study at zero temperature will be per- 
formed using exact diagonalization techniques on finite 
clusters. Our starting point is the fact that numerical 
studies of the pure t—J model on the square lattice 
have shown that a two-hole bound state is formed for 
J > J c ~ 0.3 t in the bulk limit and with d x 2_ y 2 in- 
ternal symmetry. |ll],^(J. These studies have also given 
strong unbiased indications of d x 2_ y 2 superconductivity 
at quarter filling (n = 0.5) in the vicinity of phase separa- 
tion, (2l]j2^] although they are less conclusive near half- 
filling. Our objective is also to ellucidate which is the 
pair-breaking mechanism specially for the case of non- 
magnetic impurities. 

The Hamiltonian of the t — J model in the presence of 
impurities is: j[7| 

H= - t (cjo^o- + Co/*) 

+ J ^ ( S * ' S J ~ \ n i n j) 

<i,j>,i,j^io 

+j' E ( s *° • s ^ - I (i) 

<io,j> 

where cL, is an electron creation operator at site i with 
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spin a with the constraint of no double occupancy, m = 
n$ + n$4 is the electron number operator. The impuri- 
ties are located at the sites {io}- Here t is the hopping 
parameter, which we choose as the scale of energies, J 
is the AF exchange interaction and J' is the magnetic 
coupling between the impurity spin and the electrons in 
nearest neighbor (NN) sites. The special case J' = cor- 
responds to an inert site describing a Zn impurity. ]2j ] 
We adopt periodic boundary conditions throughout. 

As a preliminary study, it is interesting to investigate 
the effect of the magnitude of the impurity spin iS$ . Fol- 
lowing previous work, we define the hole-impurity bind- 
ing energy Q , 

E B ,i = e(lh, li) - e(0h, li) - e(lh, Oi), (2) 

where e(Nh,Mi) = E(Nh, Mi) - E(0h,0i), E(Nh, Mi) 
is the ground state energy of the system with N holes and 
M impurities. For the case of a spin 1 impurity embedded 
in small clusters, bound states are found (i.e. Eb.i < 0) 
in all symmetry channels provided | J'/ J | does not ex- 
ceed small critical values. This is qualitatively similar to 
the spin 1/2 impurity as clearly seen in Fig. [l] show- 
ing a comparison of Eb,i calculated on a 20-site cluster 
in the d-wave channel for spin 1/2 and spin 1 impurities. 
Note that the coupling J' is more effective to destroy the 
bound-state for larger impurity spin. This can be simply 
understood from a simple argument relating the binding 
energy with the number of broken bonds in the AF back- 
ground |l(| (this also applies to pairing of holes |l^]). In- 
deed, while the average magnetic energy per bond for the 
bonds not connected to the impurity ]24[ ] depends weakly 
on J', the magnetic energy for the bonds connected to 
the impurity essentially scales like J' Si so that the ef- 
fective short range attractive potential for holes weakens 
for increasing J' or Si . In the following, we shall restrict 
ourself to the lowest energy sector. Moreover, since from 
Fig. the results for S = 1/2 and S = 1 are very similar 
for the region of interest J' > 0, we shall consider only 
the case of spin-1/2 impurities. 

In order to determine the binding of holes in pairs in 
the vicinity of impurities let us start examining the fol- 
lowing combinations of ground state energies: 

E B ,2= e(2h,0i) - 2e(lh,0i) 

E' B 2 = e(2h, Oi) + e(0h, li) - e(lh, li) - e(lh, Oi) (3) 
E B ,3= e{2h, li) - e(0h, li) ~ e(2h, Oi) 

The physical meaning of these quantities is quite obvi- 
ous. E B ,2 is the usual binding energy for the pure sys- 
tem. Eb.3 corresponds to compare the state where 
the two holes are trapped to the impurity with the state 
where the holes move away from the impurity. E' B 2 is the 
energy difference between the state with a hole pair away 
from the impurity and the state with one hole bound to 
the impurity. Thus, for the holes forming a bound pair 
not trapped by the impurity the following inequalities 
should hold: E B .i < 0, E' B 2 < 0, and E B 3 > 0. 



The results for E B ,i, E' b 2 and E B ,3 obtained for the 
4x4 lattice and for the case of a nonmagnetic impurity 
J' = are shown in Fig. ^(a). Very similar results are 
also obtained for the 18 sites tilted cluster. It can be seen 
that the conditions for the existence of a hole pair not 
bound to the impurity are satisfied in the range 0.2 < 
J/t < 0.5. For a better understanding of the physical 
situation that emerges from the study of E B ,2, E' b 2 and 
Eb.3, it is important to examine also the hole- impurity 
binding energy E B ,i as well as 

E BA = e(2h,2i)-2e(lh,U)- (4) 

Notice that Eb.a a s E B 2 reduces to E B ,2 in the ab- 
sence of impurities. Eba > and E B ,\ < means 
pair-breaking effect, each hole been trapped by an im- 
purity. Both quantities have been added to Fig. ^(a). 
Eb.i is negative for J > 0.2 implying that there is a 
hole-impurity bound state. It is interesting to note that 
taking into account the results for Eb.2, E' b 2 and -E'b, 3 
discussed above this hole-impurity bound state only ap- 
pears for J > 0.5 when a second hole is added. However, 
as we shall show below, the tendency to form a hole- 
impurity bound state is the main source of suppression 
of pairing. Also interesting is the result of Eba indicat- 
ing that two impurities split the hole pair for J < 0.3. 

The most important results for their connection to the 
experimental results observed in the cuprates are shown 
in Fig. §(b). In this figure, we show E' B 2 and Eb.3 ob- 
tained for the 4x4 cluster for several values of J'/ J. 
We also add E B ,2 for comparison. It can be seen that 
E' B 2 is suppressed as J' /J decreases from 1, correspond- 
ing to the pure case except for the condition of exclusion 
of holes at the impurity sites, to which corresponds to 
the nonmagnetic case. This behavior of the binding en- 
ergy in the presence of impurities is consistent with the 
experimentally observed one. Eb,3 also decreases in such 
a way that the interval of J/t where the bound state of 
holes exists and is not trapped to the impurity narrows 
down as J' /J decreases. 

Additional evidence of the strongest pair-breaker effect 
of nonmagnetic impurities in the t — J model comes from 
the study of the quasiparticle weight of pairs and pairing 
correlations with d x 2_ y 2 symmetry. The quasiparticle 
weight of pairs Z2h is defined as 

I < ^|A|^ > I 
2h (< ^lAtAI^ >)V2 [0> 

where A = -^J2i^i A». The pairing operator at site 
i 7^ io is 

A* 

where the sum extends over the NN of site i and are 
the form factors that determine the pairing symmetry. 
Note that when i corresponds to one of the four NN sites 
of the sum over /j, is restricted to only i + fi ^ iq. 
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The rest of the notation is standard. pp| , p5t The results 
obtained in the 4x4 lattice for d x 2_ y 2 symmetry as a 
function of J/t and J'/J = ,0.5 and 1.0. are shown in 
Fig. ||(a). It can be seen that the largest reduction of 
Z211 corresponds to the J'/J = case. Notice from the 
result for J'/J = 0.5 that this reduction is not linear in 
J' I J . The pairing correlations are defined as 

where, in the presence of an impurity, the sum is re- 
stricted to i and i + r =/= iq. It is well-known that close 
to half-filling the pairing correlations at large distances 
are quite suppressed, []l8|| For this reason, in Fig. [|(b) 
we show Ca(t) for r > 1 on the 4x4 lattice at quarter 
filling with J/t = 3 where they are maximum in the pure 
system, plj The tail in C\(r) is mostly suppressed for 
J'/J = 0. The results for J'/J = 0.5 indicate again a 
nonlinear behavior of pairing as a function of J'/J. 

In order to understand the origin of the pair breaking 
mechanism let us examine the quantity W p defined as the 
probability of finding the pair of holes in NN sites. In Fig. 
0(a) we show W p as a function of J/t for the pure system 
(full circles) and in the presence of a single impurity with 
J'/J = (full squares). It is also shown for various J'/J 
the quantity W p computed with the constraint that the 
holes can not occupy the NN sites of the impurity. For 
the case of J'/J — 0, by comparing the total W p with the 
restricted one it can be concluded that there is a large 
contribution to W p coming from the hole pair bound to 
the impurity. This contribution becomes smaller as J'/J 
is increased from zero. 

To complete the characterization of the effects of im- 
purities on the ground state of the t — J model we briefly 
discuss their effect on the magnetic order |l5|] and on the 
dynamical magnetic fluctuations. The dynamical spin 
structure factor at the AF wavevector (71", w) shown in 
Fig. [|(b) exhibits, in comparison to the pure case jTsLEfl, 
new structures at low energies. Spectral weight is trans- 
fered from the peak at w ~ 0.8J, characteristic of the 
pure two-hole doped 4x4 system, to lower energies of 
order J'. This new resonance can be interpreted as the 
singlet-triplet excitation energy of the singlet impurity- 
hole bound-state. These results can be related to recent 
neutron scattering experiments which have reported that 
the magnetic pseudogap at the AF wavevector disappears 
with Zn doping. |^7|l 

In summary, we have studied using exact diagonaliza- 
tion of the 2D t — J model in the presence of magnetic 
(S = 1/2) and nonmagnetic (S — 0) impurities. We have 
shown through an analysis of binding energies, pairing 
correlations and dynamical pair susceptibility that the 
nonmagnetic impurity has a stronger pair breaker effect 
than the magnetic impurity. These results arc in agree- 
ment with experimental findings in Zn (S — 0) and Ni 
(5 = 1) doped cuprates. This agreement reafirms the 
ability of microscopic strongly correlated electronic mod- 



els to describe the physics of CuC>2 planes in the cuprates. 
We have also given some indications of the possible es- 
sential role of the binding of holes to impurities in sup- 
pressing superconductivity in these systems. 
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comments. Part of the numerical calculations were done 
at IDRIS (France). 
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FIG. 1. Eb,i vs J'/J for d-wave orbital symmetry, for a 
spin 1/2 impurity (solid squares) and a spin 1 impurity (open 
squares) . 

FIG. 2. a) Eb,2 (open squares), E' B 2 (open diamonds), 
Eb,3 (solid circles), Eb.\ (solid squares) and Eb,a (solid trian- 
gles), defined in the text, as a function of J/t and for J' — 0. 
b) Eb,3 (open symbols) and E' B 2 (solid symbols) vs. J/t. 
The results correspond to J' = (squares), J' = J/3 (cir- 
cles), J' = 2 J/3 (triangles) and J' = J (diamonds). We also 
show for comparison the result of Eb,2 (crosses). 

FIG. 3. a) The quasiparticle spectral weight Z 2 h as a 
function of J/t for different values of J'/J. b) The pairing 
correlation Ca{t) at quarter filling for different values of J'/J. 
The results for the model with no impurities are shown with 
solid circles. 

FIG. 4. a) The probability W p vs. J/t. The results for 
the pure system are shown with solid circles while with solid 
squares we show the results with one impurity in the system 
and J' = 0. The results for W p with the exclusion of the 
NN sites of the impurity are shown for various values of J'/J 
with open symbols as indicated on the plot, b) Dynamical 
spin structure factor (arbitrary units) in the two hole doped 
system for J — 1.0, and J' = (solid line), 0.25 (dashed line), 
0.5 (dotted line) and 1.0 (thick solid line). 
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